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The random antiferromagnetic two-leg and zigzag spin-1/2 ladders are investigated using the 
real space renormalization group scheme and their complete phase diagrams are determined. We 
demonstrate that the first system belongs to the same universality class of the dimerized random 
spin-1/2 chain. The zigzag ladder, on the other hand, is in a random singlet phase at weak frustration 
and disorder. Otherwise, we give additional evidence that it belongs to the universality class of the 
random antiferromagnetic and ferromagnetic quantum spin chains, although the universal fixed 
point found in the latter system is never realized. We find, however, a new universal fixed point at 
intermediate disorder. 
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INTRODUCTION 



One dimensional spin systems have been extensively 
studied over the last several years and a fairly deep un- 
derstanding of their possible phases and corresponding 
physical behavior has emerged. Although actual real- 
izations are restricted to systems with a high degree of 
anisotropy, one of the main motivations for these studies 
is the possibility of cataloging their universality classes. 
This is specially tempting in the case of disordered sys- 
tems, since specially powerful methods can be employed 
in one dimension, which are able to expose not only av- 
erage values but full distribution functions. One of the 
most studied random one-dimensional spin systems is the 
random antiferromagnetic (AF) spin-1/2 chainii*^ Mak- 
ing use of the real space renormalization group (RSRG) 
method of Ma, Dasgupta and Hup* it has been shown 
that, for any amount of uncorrelated disorder^ the low- 
energy physics of this system is governed by an infi- 
nite randomness fixed point^ The approach to this fixed 
point is characterized by the formation at decreasing en- 
ergy scales of random singlets between widely separated 
spins. In this random singlet (RS) phase several phys- 
ical properties are universal and knowni^ For example, 
the spin susceptibility \ ~ l/Tlog 2 T, and the spin- 
spin correlation function Cy = (S,; • Sj) is such that its 
mean value Cy ~ (— 1)* J / |i — j\ 2 , while the typical one 

\ c a\typ ~ e M-V\i-T\)' 

Other random systems have also been analyzed by 
these methods, among which two are of special inter- 
est to us. One is the random dimerized AF chain with 
different distributions for odd and even links 4 The pres- 
ence of a gap in the clean version of this system provides 
protection against the introduction of disorder. There- 
fore, for weak disorder, the system retains a gap (or a 
pseudogap) and a spin susceptibility which decreases to 
zero with decreasing temperature. For strong enough 
disorder, however, the pseudogap is destroyed and the 
susceptibility diverges as a power law. The power law 
exponent is non-universal and varies continuously with 
the disorder strength, characterizing a Griffiths phase. 



The other system of interest is the random chain with 
both AF and ferromagnetic (FM) interactions*^ It has 
been shown, using a generalization of the RSRG pro- 
cedure, that the low-energy behavior of this system is 
governed by the formation of long clusters with large to- 
tal spins. For weak disorder, this large spin (LS) phase 
is universal. The magnetic susceptibility diverges like a 
Curie law, and the specific heat vanishes like T 1 I Z \ lnT\, 
with 1/z « 0.44. However, for strong disorder, the LS 
phase is no longer universal: although the magnetic sus- 
ceptibility is still Curie-like, the specific heat exponent z 
varies continuously with disorder. 

More recently, following the discovery of spin ladder 
materials, 7 attention has been drawn to coupled pairs 
of spin chains. The random two-leg and zigzag ladders 
are the most studied of them^^IS The two-leg ladder 
is known to possess two phasesi^ a gapped phase, with 
vanishing spin susceptibility x (T), and a Griffiths phase, 
where \ (T) diverges as a power law with continuously 
varying exponents as T — » 0. The zigzag ladder, on the 
other hand, is topologically equivalent to single chains 
with both nearest neighbor (nn) and next nearest neigh- 
bor (nnn) interactions. Initial studies were confined to 
fairly small systems and concluded in favor of the exis- 
tence of only one phase of the Griffiths type£ Subsequent 
investigations identified the presence of a small RS phase 
for weak nnn interactions*^ Furthermore, the formation 
of large effective spins, with FM and AF nn and vanish- 
ing nnn couplings at the late stages of renormalization 
were taken as an indication of a LS phased 

The main purpose of this work is to establish in detail 
the effective low-energy equivalence between the random 
two-leg and zigzag ladders with the three random spin 
chains mentioned above. We do so by means of the RSRG 
method. In particular, we will show that at low energies 
the random two-leg ladder is equivalent to the dimer- 
ized AF spin-1/2 chain. In the course of this analysis, 
we will construct the full phase diagram of this model 
and confirm the existence of only two possible phases. 
Moreover, we will show that there are indeed two possi- 
ble low-energy behaviors in the zigzag ladder. For small 
and weakly disordered nnn couplings, it is equivalent to a 
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random AF spin-1/2 chain with its RS phase. If the nnn 
interactions are stronger, however, the system presents 
all the characteristic features of the LS phase of a ran- 
dom chain with both AF and FM couplings. We thus 
confirm the results of Ref. ^] by showing the same scal- 
ing of the effective cluster sizes and total spin values with 
energy as seen in those systems. However, the zigzag lad- 
ders can never exhibit the universal behavior observed in 
those systems at weak disorder. This is demonstrated 
by an exhaustive investigation of the dependence of the 
dynamical exponent z on the shape and strength of the 
disorder distribution. Nevertheless, we do find different 
universal regions where z remains fixed while the disorder 
is varied. The totality of our results thus serves to show 
that the universality classes of random spin ladders can 
already be found in simpler systems with nn interactions 
only. 

The paper is organized as follows. In Sec. [H] we in- 
troduce the models we studied. The numerical results 
of the RSRG procedure of both the conventional two-leg 
and zigzag ladders are presented in Sec. IIIII The unex- 
pected behavior of the dynamical exponent of the zigzag 
ladders is explained in Sec. IIVI We end with some final 
conclusions in Sec. [3 Appendix El discusses our method 
of calculation of the dynamical exponent z and a brief 
discussion of the case of correlated disorder is given in 
Appendix IbI 



In order to study these systems we employed the RSRG 
method introduced by Ma, Dasgupta and Hui Its dec- 
imation steps consist in isolating the strongest bond of 
the system (0), keeping only the lowest energy level of 
the bond, and renormalizing the remaining interactions 
by perturbation theory. The new renormalized coupling 
constants can be either ferromagnetic or antiferromag- 
netic. The details of the procedure have been extensively 
discussed in the published literature^Siifi and will be 
skipped here. As the largest energy scale is lowered from 
its initial value Qq — max(J ,ifo), an effective distribu- 
tion of coupling constants is generated, which eventually 
flows towards a fixed point distribution. The low energy 
behavior of the system is governed by the remaining 'ac- 
tive' non-decimated spins at the scale of interest ! 1 ' 2 " 1 ? 

When there are initially only nn interactions, further 
neighbor interactions are never generated ji* 2 ^ By con- 
trast, in our case, the range of effective couplings rapidly 
increases as the RSRG is iterated (see Section|^for more 
details). However, as the fixed point is approached, in- 
teractions beyond nearest neighbors become extremely 
weak. If we then neglect interactions weaker than a 
certain upper bound (O m j„ w 10 _ 200 f2 )±i the effective 
range actually extends only as far as the nearest neigh- 
bors. Thus, at the final stages of the RSRG, the ladders 
renormalize to effective nearest neighbor spin chains. 



II. THE MODEL 

We consider the following Hamiltonian, which de- 
scribes an AF spin-1/2 chain with nearest neighbor (nn) 
and next nearest neighbor (nnn) interactions 



H 



N-l 

E 



JV-2 



K& • S 



i+2) 



(1) 



where Si is a spin-1/2 operator, N is the total number of 
spins, and Ji > e Ki > are the nn and nnn random 
coupling constants, respectively. If Ji = for even i, this 
is the two-leg ladder Hamiltonian; if Ji is in general non- 
zero, we have the zigzag ladder. The non-zero coupling 
constants Ji and Ki are in general independent random 
variables (see, however, Appendix [Bj ■ We take them to 
be respectively distributed in a power-law fashion (unless 
otherwise noted): 



Pj{J) 



1 — a / Jo 

~TT It 



< J < Jo, (2a) 



Pk{K) = L^f^) . „. /v K„. (21,) 



The exponent a(0<a<l)isa measure of the dis- 
order strength and the ratio of cutoffs A = K Q /J gives 
the typical relative strength between nnn and nn inter- 
actions. 



III. NUMERICAL RESULTS 

In this section we show the numerical results obtained 
from the iteration of the RSRG for the two-leg and the 
zigzag ladders. 



A. Two-leg ladders 

We first focus on the two- leg ladders {Jn = 0). In our 
simulations we used chain lengths up to N — 200,000. 

In Fig.^ we show the behavior of the fractions of nn, 
nnn, and 3rd nn bonds as functions of the energy scale 
O. In addition, this Figure also shows the fraction of 
'active' spins greater than 1/2 and the fraction of anti- 
ferromagnetic coupling constants as functions of fl. The 
first thing to note is that the only significant couplings at 
the lowest energy scales (< 10 -3 5 Qq) are nn couplings. 
Furthermore, these remaining couplings are all antiferro- 
magnetic. This has been verified for all values of a and 
A. 

Another important feature of the approach to the fixed 
point is the difference between the distributions of the 
odd bonds and the even ones. We have checked that, at 
the lowest energy scales, all the even bonds vanish while 
the odd ones are distributed according to 



Pjodd(J) ~ J- W/Z , 



(3) 



where z is the dynamical exponenti&i^ Thus, the sys- 
tem renormalizes to a collection of free effective random 
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Figure 1: The behavior of the fractions of nearest-neighbor, 
next-nearest-neighbor, and third-nearest-neighbor coupling 
constants, the fraction of spins greater than 1/2, and the frac- 
tion of antiferromagnetic couplings as functions of the energy 
scale fi. The initial number of spins is N = 200, 000, the dis- 
order strength is a = 0.3, and the ratio A = 1.0. The RSRG 
is iterated until only nearest-neighbor couplings are left. The 
data are averaged over 100 samples and the relative error is 
less than 2%. 



dimers. The average magnetic susceptibility is then given 
byM 

If z < 1, the density of states is suppressed at low energies 
and the system has a 'soft gap' (or a pseudo-gap). This is 
a remnant of the Haldane-type gapA^ of the clean two-leg 
ladder We therefore call it a disordered Haldane phase. 
Otherwise, the system shows no such suppression and is 
in a gapless (Griffiths) phased We have determined the 
value of z, for different disorder strengths (see the Ap- 
pendix El f° r details on the method of computation of 
z). As shown in Fig. El for A = 1, a transition between 
these two phases occurs when a ~ 0.3. In both phases 
the system is strongly dimerized. Within the length scale 
of the size of the effective dimers, the correlations decay 
as a power law. We believe this decay is similar to the 
case of random Heisenberg chains (~ l/r 2 ), since there 
are equal contributions coming from even and odd bond 
decimations. At larger length scales, the correlations are 
strongly (exponentially) suppressed. This overall behav- 
ior is apparent in Fig. 5 of Ref.0. The crossover between 
the two regimes is governed by the effective dimer sizes, 
which are primarily determined by A not a. 

Calculating the dynamical exponent z for various val- 
ues of A and a, one can construct the phase diagram of 
the two-leg ladder (Fig. EJ) • In the transition line, the dy- 
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Figure 2: Variation of the dynamical exponent z with the 
disorder strength a for A = 1.0. N = 200,000 is the initial 
number of spins in the ladder and the RSRG is iterated un- 
til there are only nearest-neighbor bonds left. The data are 
average over 5 samples and the relative error is about 3%, 
comparable to the symbol size. 



namical exponent equals one and the low-energy density 
of states and the magnetic susceptibility as T — > are 
both constant. 

From the inset of Fig. we note that for A < 10~ 2 
the system is always in the Griffiths phase. We can un- 
derstand this in the following manner. In the limit of 
A — > 0, the system reduces to a collection of discon- 
nected dimers (the 'rungs' of the ladder), whose cou- 
plings are distributed according to Eq. (l2al) . Therefore, 
z = (1 — a) 1 > 1 and the system is always in a Grif- 
fiths phase. For small A, this behavior is preserved. For 
10~ 2 < A < 1.53, the Haldane-type gap of the clean 
system gives rise to a 'soft gap' upon the introduction 
of weak disorder. For large disorder, the Griffiths phase 
re-emerges, with a diverging non-universal magnetic sus- 
ceptibility. For A > 1.53, only the Griffiths phase exists. 
This behavior smoothly connects with the A — » oo limit 
of two disconnected random Heisenberg chains, which is 
governed by the infinite randomness fixed point^ For- 
mally, this limit corresponds to z — > oo. We note that 
for any value of A, there are only two types of phases. 8 "® 

We stress the high degree of similarity between the dis- 
ordered two-leg ladder and the dimerized random antifer- 
romagnetic chain, 4 ' which also shows analogous phases to 
the ones discussed here. Indeed, in both cases the RSRG 
flow leads to a fixed point with nn antiferromagnetic in- 
teractions only but with different distributions for even 
and odd bonds. Thus, the two systems clearly belong to 
the same universality class. 
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Figure 3: Phase diagram of the random two-leg ladder. The 
disordered Haldane phase is characterized by a 'soft gap' (or 
pseudo-gap), whereas in the Griffiths phase, the disorder com- 
pletely destroys this remnant of the Haldane gap. In both 
phases the system is dimerized. The data error is about the 
size of the symbols. 



B. Zigzag ladders 



Figure 4: The fraction of nearest-neighbor, next-nearest- 
neighbor, third-nearest-neighbor, the fraction of spins greater 
than 1/2, and the fraction of AF couplings for zigzag ladders 
as functions of the energy scale Q. The initial number of spins 
is iV = 640, 000, the disorder strength is a = and the ratio 
A = 1.0. The renormalization group is iterated until all the 
bonds in the system are between nearest-neighbors only. The 
data are averaged over 5 samples. For ft/Qo > 10 -10 the data 
relative error is less than 2%, otherwise, it is less than 8%. 



We now consider the zigzag ladders, with coupling con- 
stants distributed according to Eq. (2J. We studied lad- 
ders with initial lengths of N = 640, 000. 

Yusuf and Yang 1 ? have calculated the ratio of the av- 
erage nn interactions to that of further neighbor inter- 
actions in zigzag ladders with correlated disorder. They 
have shown that this ratio increases tremendously as the 
energy scale is lowered. Indeed, if we plot the fraction of 
nn, nnn, and 3rd nn interactions as functions of fl, for 
a = and A = 1 (see Fig. EJl , we can see that only nn 
bonds survive at the lowest energies (« 10~ 15 f2o)- This 
is similar to the case of two- leg ladders, as shown in Sec- 
tion |^EJ and was verified in our simulations, irrespec- 
tive of the values of a and A. However, in contrast to the 
two-leg ladder case, the nn bonds are not all antiferro- 
magnetic. The asymptotic number of ferromagnetic and 
antiferromagnetic bonds is about the same, as shown also 
in Fig. 31 Finally, as also pointed out in Ref. 13, there 
is a rapid proliferation of spins greater than 1/2. This 
is typical of the so-called Large Spin (LS) phase found 
in disordered Heisenberg chains with both FM and AF 
interactions^ As we will see, this similarity is not fortu- 
itous. 

One of the distinguishing features of the disordered 
chains with both FM and AF interactions is the re- 
lation between the average spin size and the average 
cluster length. The average cluster size is related to 
the energy scale through the dynamical exponent z by 



< I >~ p~ x ~ _1 / z and the average spin size grows 
with the lowering energy scale according to < S >~ Q~ K . 
Both behaviors are observed in the zigzag ladders as 
shown in Fig. As first pointed out by Westerberg et 
al.fi the two exponents are related by l/z = 2k in random 
chains with both FM and AF interactions. This follows 
from the fact that the main decimation process is not 
singlet formation but rather the formation of large spins 
from the random addition and subtraction of spin pairs. 
As a consequence, the cluster growth is characterized by 
a random walk in 'spin size space'. Again, we find this 
relation also holds for the zigzag ladders (k « 0.069 and 
l/z « 0.14 for the chain of Fig.^J. This is evidence that 
the zigzag ladders belong to the same universality class 
of the disordered chains with FM and AF interactions. 

Our simulations also show that the fixed point distri- 
butions of FM and AF link excitation gaps are the same, 
with a characteristic power law dependence^ 

P AF (A) = P F (A) = P (A) ~ A- 1+1 / z , (4) 

where z is the same dynamical exponent obtained from 
the relation between < I > and fl (see below for an ex- 
planation) . 

An interesting feature of the zigzag ladders is the fol- 
lowing. In general, the total fixed point distribution of 
link excitation gaps is a linear combination of AF and 
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Figure 5: The average cluster size < I > and the average spin 
size < S > as functions of the energy scale Q. Prom a direct 
fit we get < I >~ fi-°- 14 and < S >~ ft" - 069 . The initial 
length of the system is N = 640, 000, a = and A = 1.0. The 
renormalization group is iterated until only nearest-neighbor 
couplings are left. The data were averaged over 5 samples and 
the sample to sample variation is less than 5%. 



FM contributions, namely. 
P(A) = ' 




I-I/ZAF 



1 — 1/ ZFM 



Z F M \ A 



where the first and second terms come from the AF and 
FM link distributions, respectively, and x is the fraction 
of AF couplings. In our simulations, the fraction of AF 
bonds is x « 0.53 in the LS phase. The relation between 
length scale and the energy scale is given by£ 



dp 

dn 



= P(A = fi)p. 



(5) 



The coefficient in front of P (A = f2) is taken to be 1 
because the main decimation process replaces 2 spins by 
1. Solving Eq. (jj we get 



P ' 



n 1 ^ 



with 



l/z = x/z AF + (1 - x)/z FM - 



(6) 



This general relation seems to have been unnoticed in 
previous studies. In particular, we note that weakly dis- 
ordered chains with both AF and FM couplings have 
~ 0.56, 1/zaf ~ 0.30, and x ~ 0.63^ so that 



l/z w 0.40 which is reasonably close to the value of 
0.44 ± 0.02 found directly in the simulations of Ref. 0. 
In contrast, in the zigzag ladders, we have found Zaf = 
zfm, irrespective of the value of x (x 0.53 in the LS 
phase). As a result, zaf — zfm = z. 

Now, we would like to explore the disorder dependence 
of the dynamical exponent z. Fig. illustrates the be- 
havior of z as a function of the disorder parameter a 
for various values of A. Note that the smallest value 
of z is 1/0.15. Westerberg et al. have found that for 
a range of weak disorder strengths, random chains with 
both FM and AF couplings are characterized by a uni- 
versal value of zu = 1/0. 44^ Strongly singular disor- 
der distributions, however, show non-universal disorder- 
dependent values of z. As can be seen from Fig. for 
initial power-law distributions of couplings, zigzag lad- 
ders always have z > zjj and are never in the basin of 
attraction of the universal behavior. 

Zigzag Ladder 
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Figure 6: Variation of the dynamical exponent z with the 
disorder parameter a for different values of A. From top to 
bottom, the solid lines with circles correspond to A = 1.0, 0.6, 
0.4, 0.2, 0.1, 0.05, and 0.0025, whereas the dashed lines with 
stars are for A = 0.8" 1 , 0.6~\ 0.4" 1 , 0.2" 1 , 0.1 -1 , 10 3 , and 
10 6 . = 640, 000 is the initial number of spins in the ladder 
and the renormalization group is iterated until the bonds in 
the system are nearest-neighbor only. The data relative error 
is less than 10%, about the size of the symbols. 

For a < 0.6 the behavior of l/z is linear in a 
l/z = aa + b. 

Note that, for A > 1.0, a = and for 1.0 < A < 1/0.4 
the system flows to a new universal fixed point where 
l/z — 0.15 ±0.02. This is similar to the random AF and 
FM chains but the value of the dynamical exponents are 
different. We note from the general trend of the curves 
in Fig. that the origin of this universal behavior seems 
to be the non-monotonic behavior of z as one goes from 
A < 1 to A > 1 (see also Fig.dbelow). This behavior, on 
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the other hand, can be understood from simple physical 
arguments as shown below in Section Hvl 

For A < 0.4 all the lines have the same slope a = 
0.039 ± 0.001 and the intercept b varies in a logarithm 
manner with A, i.e., b = b\ + 62 In A, where b\ = 
0.12 ± 0.01 and 6 2 = 0.020 ± 0.001. Thus, there is a 
phase transition at the value of a where the dynamical 
exponent diverges, a — a c (A) = — (62 hi A + 61) /a. For 
a < a c (A) , the system is governed by an infinite ran- 
domness fixed point where the magnitude of the spins 
does not grow and the FM couplings vanish. Indeed, the 
presence of such a RS phase was observed by Yusuf and 
Yangitt by direct calculation of z. Hence, we conclude 
that the system presents two phases. In the greater part 
of phase diagram it is in the LS phase, with a new univer- 
sal fixed point in the region 1.0 < A < 1/0.4 and a < 0.6. 
In addition, there is a tiny region where the system is in 
a RS phase. 

For a > 0.6, all the curves converge to the point 
a = 0.95 and z — 1/0.05. Note in this connection that 
a = 0.95 corresponds to an initial distribution which is 
the same as the fixed point one if z = 1/0.05, cf. Eqs. 12;ill 
and P. This is a reflection of the inability of the deci- 
mation procedure, which is dominated by first order per- 
turbation theory steps, to generate distributions which 
are more singular than the initial one. 
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distribution with a gap', where the initial distributions 
are uniform and have support in Jq — S < J < Jq and 
K - 5 < K < K (we take max{J , K Q } = 1). We 
have explored several values of S and found z to be al- 
ways greater than zjj . This is seen, for example, in the 
extreme case of 6 = 0.05 in Fig. Q where we plot 1/z 
as a function of A. Note that there is a phase transi- 
tion at A = 0.4. To the left of this point, the system 
is in a RS spin-1/2 phase and the dynamical exponent 
diverges z ~ — InO. For A > 0.4, the system is in the 
LS phase. The smallest value of the dynamical exponent 
is z = 1/0.22 (still greater than the universal z\j found 
by Westerberg et. al) and for A > 10, z saturates at 
z = 1/0.11. We conclude that the LS phase of the zigzag 
ladders is in the same universality class of random FM 
and AF chains with strong disorder but it does not reach 
the universal region found in those chains when the dis- 
order is weak. It should be remembered that the clean 
system is gapless if A < 0.24 but spontaneously dimer- 
ized and gapful if A > 0.24ii& The topological nature of 
the dimerized state, however, makes it unstable with re- 
spect to any finite amount of disorder and the random 
system is always gapless pii 
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Figure 7: Variation of the dynamical exponent z with the 
ratio A for initial 'box distributions with a gap' (5 = 0.05, see 
text for details). N = 200,000 is the initial number of spins 
in the ladder and the renormalization group is iterated until 
the bonds in the system are nearest-neighbor only. The data 
relative error is less than 5%, about the symbol size. 

One might think that the universal behavior of disor- 
dered chains with AF and FM interactions where z = zy 
could be realized in zigzag ladders with disorder distri- 
butions which are not as singular as Eq. 1J2J. This is not 
the case, however, as exemplified by the case of a 'box 



Figure 8: Phase diagram of the random zigzag ladder. The 
Large Spin phase dominates most of the parameter space. 
There is a universal region (the shaded area) for 1 < A < 2.5 
and < a < 0.6 characterized by z ~ 1/0.15. There is also 
a tiny region for very small values of A (see the inset), where 
the system is in a Random Singlet phase. 

Based on the behavior of the dynamical exponent z, we 
can determine the phase diagram of the random zigzag 
ladder, as shown the Fig. 03 The Large Spin phase is 
dominant in most of the (a, A) parameter space. In this 
phase, the low energy physics of the system is governed 
by a fixed point where the mean size of the spin clus- 
ters grows when the energy scale is lowered, and they 
are weakly coupled. The distribution of link excitations 
gaps is not universal except in the shaded region. This 
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universal region in the Large Spin phase is character- 
ized by a dynamical exponent z « 1/0.15. The ther- 
modynamic properties are well known: 5 the specific heat 
C ~ T x l z \ InT and there is a Curie-like magnetic suscep- 
tibility x ~ T _1 . The average spin-spin correlation func- 
tion decays in a logarithmic manner^ In addition to this 
phase, there is a tiny region where the system is in a RS 
phase (where z — > oo, see the inset of the Fig. EJ. Here, 
the low energy physics is governed by a universal infinite 
randomness fixed point. The specific heat C ~ —1/ In 3 T, 
the magnetic susceptibility \ ~ 1/ (T In 2 T) , the mean 



couplings are 



i(r) 



', and the typical corre- 



exp | — r 1//2 }i^ This phase was 



correlation function C, 
lation function C t (r) r 

previously identified in Ref . Oil, although its precise loca- 
tion differs somewhat from our results. We attribute this 
difference to finite-size effects and the different methods 
used for the characterization of RS behavior. 



IV. ANOMALOUS DISORDER DEPENDENCE 
OF z IN THE ZIGZAG LADDERS 

An intriguing aspect of the data shown in Fig. U3 is 
the fact that z decreases as the initial disorder strength 
a is increased, for A < 1 and a < 0.6 (solid lines of 
Figure EJ- This behavior is unexpected since the weaker 
the initial disorder is, the stronger is the final effective 
one. As we will show, this anomalous behavior can be 
simply understood by analyzing the limiting behavior as 
A^O. 
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Figure 9: Schematic decimations at the earliest stages of the 
RG. 



Let us start by writing down the effective couplings 
generated by the decimation procedure in its earliest 
stages when all the spins are spin- 1/2 and the couplings 
are AF (see Fig. EJ. Let us assume that the coupling 
between spins £3 and S4 is the largest one in the sys- 
tem, i.e., SI = J34. After the decimation step, both spins 
are removed from the system because they are effectively 
frozen in a singlet state. The remaining spins are cou- 
pled with renormalized interactions as follows. The nn 



J12 

J25 
J56 



= J12- 



= J; 



if 13 ( J23 ~ ^ 24) 
20 

( J 23 - K 24 ) ( J 45 - K 35 ) 
20 

i^46 (^45 — K 35 ) 



56 



20 



the nnn ones are 



Kib = 
K 26 = 



if 13 (J45 ~ if 35) 
20 

if 46 (J23 ~ if 24) 
2J1 



and a third-nearest-neighbor coupling between spins Si 
and Sq is generated 



'16 



ifi3if 46 
20 



First, we analyze the limit A < 1. In this case, it 
is easy to see that the renormalized nn, nnn, and 3rd 
nn couplings are of order 0(A°), ©(A 1 ), and 0(A 2 ), re- 
spectively. Thus, if we neglect the 3rd nn coupling, the 
original form of the Hamiltonian is recovered, and the 
nn couplings remain stronger than the nnn ones. No FM 
coupling is generated and as the energy scale is lowered 
the nnn couplings vanish faster than the nn ones. The 
system flows to a random AF spin- 1/2 chain in a RS 
phase. This scheme breaks down when a nn FM cou- 
pling appears leading to frustration. Among the three 
nn renormalized couplings, J25 is the most likely to be 
FM with probability given by 



P Jm < 



2P ( J 23 < K 2i ) P ( J45 > if; 



35; 




J Pk (if) dK 
jl-^ Pj{J) Px (if) dK 



dJ} x 



dJ 



A 1 -" f l - ^A 1 - ! 



As expected, this probability is greater, the greater A 
is, since A governs the strength of frustration. But note 
also that it increases with a. As more FM couplings 
are produced, the RG flow tends to deviate from the RS 
phase. As a result, the dynamical exponent will tend 
to be smaller. This explains the anomalous behavior we 
have found. 

In the opposite limit (A ^> 1), this picture no longer 
holds. Suppose that O = ^35. In this case, the renor- 
malized coupling between spins S 2 and S4 is now a nn 
one 



■h 



if 2 



J23U 



45 



J: 



34; 



20 
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while the nnn coupling between Si and S4 is 

r > -^13(^45 ^ ^34) 

Kli = m ■ 

It is clear that in this case the nn coupling tends to in- 
crease and the nnn one has a better chance of becoming 
a FM coupling in the earliest stages of the RG flow. This 
will rapidly enhance the frustration of the system and 
drive it away from the RS phase. In this limit, we expect 
the RS phase only at A -1 = 0. 

The competing tendencies at small and large A lead 
to a minimum value of z (Figs. and UJ. The presence 
of this minimum partially explains the universal region 
where z is approximately constant at « 6.7 (Fig.EJ. 

A similar type of reasoning can be used to analyze the 
case of correlated disorder considered in Ref. 0. This is 
discussed in Appendix El 



the latter system has a universal fixed point at weak dis- 
order (with z w 2.27)& which is inaccessible to the former. 
The zigzag ladder, on the other hand, has a region of pa- 
rameter space with a new universal behavior: for initial 
distributions less singular than Pj (x) ~ Pk (x) ~ x~ ac , 
with a c w 0.6, and for 1.0 < A < 1/0.4, the dynam- 
ical exponent is always z ~ 6.7. This critical disorder 
strength a c is comparable to the one found in the ran- 
dom AF and FM chain^ although we see no other sim- 
ilarities between these two regimes. We have been able 
to find a rough explanation for this universal behavior 
by ascribing it to competing tendencies which create a 
shallow valley where z is approximately constant. It is 
possible that the universal regime of Ref. has a similar 
origin. 
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V. CONCLUSIONS 

We have studied random AF two-leg and zigzag spin- 
1/2 ladders using the real space renormalization group 
method of Ma, Dasgupta and Hu. We have determined 
the complete phase diagrams of these two general models 
in great detail by calculating the dynamical exponent z 
and the range of low-energy effective interactions. The 
two-leg ladders show a gapped disordered Haldane re- 
gion with a random dimer nature and a gapless Griffiths 
phase, but no random singlet phase. The zigzag ladder, 
on the other hand, can be either in a random singlet 
phase or a large spin phase. 

One of our main findings is that throughout their phase 
diagrams these two models lead to effective low-energy 
models with nearest neighbor interactions only. This sim- 
plification had been noticed before in connection with the 
zigzag ladders^ Our calculations confirm those results 
and generalize them to the two-leg ladders. Furthermore, 
by analyzing the structure of odd and even links in the 
case of the two-leg ladders and the scaling with energy 
of the large clusters of spins in the case of the zigzag 
ladders, we have been able to show the low-energy equiv- 
alence of these systems to the random dimerized AF spin 
chain and the random chain with AF and FM interac- 
tions, respectively. The low-energy equivalence between 
random systems with further neighbor yet short-ranged 
interactions to systems with nearest neighbor interac- 
tions only is likely to hold in general. However, there re- 
mains the possibility that one-dimensional models with 
longer ranged interactions are not decimated down to 
nearest neighbor spin chains in the process of the RG?i& 
It would be interesting to determine the critical range 
beyond which such simplification does not occur. 

While the zigzag ladder could be related to the random 
chain with both AF and FM chains, we have shown that 
their phase diagrams do not completely overlap. Indeed, 
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Appendix A: THE CALCULATION OF THE 
DYNAMICAL EXPONENT z 

In this appendix, we would like to comment on our 
method of calculation of the dynamical exponent z. A 
common procedure, frequently used in the literature, is 
to decimate the system until one pair of spins is left. 
One can then calculate the excitation gap of this last 
dimer (the first gap Ai). By repeating this procedure 
for different realizations of disorder one can obtain the 
distribution of first gaps. The dynamical exponent is 
obtained by fitting the behavior at small Ai to the power 
law in Eq. Q. s ' lx This is shown in Fig.EJa). Another 
option is to calculate the density of 'active' spins p at the 
scale fl. The dynamical exponent z is then obtained from 
its definition, through the relation between the length 
scale I ~ and the energy scale fl 

p-tt 1 '*. (AI) 

This second method is exemplified in Fig. iroT hl. and 
the agreement is excellent. While the methods are equiv- 
alent it should be pointed out that the second method is 
computationally much faster, as a single realization of a 
large chain is sufficient for the determination of z. In our 
calculations, we have chosen to use the second method. 
This was actually crucial in the case of the zigzag ladders, 
where the asymptotic behavior can only be obtained with 
very large system sizes rendering the first method very 
inefficient. This need for large sizes is due to the presence 
of long-range correlations in these systems^ as opposed 
to the exponential dependence of the two-leg ladders. 9 
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Figure 10: The calculation of the dynamical exponent z. (a) The distribution of first gaps Ai using 50,000 realizations for 
Y=100, 150, and 200, is fitted to log 10 P(— log 10 Ai) = const. + (l/z)logAi, and (b) z is directly calculated from log 10 p ~ 
const. + (1/z) log 10 ft, where we used Y=200,000. The data are from the decimation of a single chain. We verified that the 
variation with respect to other realizations of disorder is less than 1%. 



Appendix B: ZIGZAG LADDERS WITH 
CORRELATED DISORDER 



A', 



A 



/ J12J: 



25 



Yusuf and Yang have also considered the case of zigzag 
ladders with correlations between the random nn and 
nnn couplings 3 They have considered nn couplings dis- 
tributed according to Eq. 12al and nnn couplings given 

by 



Ki = A 



Their analysis identified a RS phase for A < 0.5, and a 
LS phase for 0.5 < A < 1.0. 

We can qualitatively analyze this case along the same 
lines as in Sec. lIVI Suppose Q — J34 in Fig. El Using the 
above definition of Ki in the earliest stages of the RG 
flow (when Q w Qq), we have 



J25 = (1 - A) 



2 J23J45 



K 26 = A' 



50 



where A' = A/(l — A). Neglecting the 3rd nn coupling 
and the renormalization of J12 and J56, the new Hamil- 
tonian has the same form as the original one, except 
for the fact that the anisotropy parameter for K15 and 
A26 has been renormalized. In the first RG decimations, 
no FM couplings can arise. However, in a second run 
there is a finite probability for FM couplings to appear 
if A' > 1. This happens if A > 0.5. The presence of FM 
couplings introduces frustration which is the mechanism 
that drives the system away from the RS phase towards 
the LS phase. 
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